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ABSTRACT 

It is possible to reliably identify white dwarfs (WDs) without recourse to spectra, 
instead using photometric and astrometric measurements to distinguish them from 
Main Sequence stars and quasars. WDs' colours can also be used to infer their intrinsic 
properties (effective temperature, surface gravity, etc.), but the results obtained must 
be interpreted with care. The difficulties stem from the existence of a solid angle 
degeneracy, as revealed by a full exploration of the likelihood, although this can be 
masked if a simple best-fit approach is used. Conversely, this degeneracy can be broken 
if a Bayesian approach is adopted, as it is then possible to utilise the prior information 
on the surface gravities of WDs implied by spectroscopic fitting. The benefits of such an 
approach are particularly strong when applied to outliers, such as the candidate halo 
and ultra-cool WDs identified by Vidrih et al. (2007). A reanalysis of these samples 
confirms their results for the latter sample but suggests that that most of the halo 
candidates are thick disk WDs in the tails of the photometric noise distribution. 
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1 INTRODUCTION 



White dwarfs (WDs) are sufficiently long-live d that they can 
be u s ed to measure the Galaxy's age (e.g., iFontaine et al.l 
I2OOII . I2OO5I and references therein) and star-formation his- 



ical probes 


'e.g.. iBergeron et al.lll997l: 


Leggett et al.l 


199S 


Harris et al. 


igggT Hodgkin et al.l 200C 


: Gates et al.l 


2004 


Harris et al. 


20081). WDs in the Galaxy's halo have also 



been invoked as potentia l dark matter candidates (e.g., 
lOppenheimer et al.ll2001al lbh. althoug h these clai ms are the 
subject of considerable debate (e.g., iReidI I2OO5I and refer- 
ences therein). 

All such investigations rely on being able to constrain a 
WD's properties: its distance, D, and its intrinsic param- 
eters (effective temperature. Toff; surface gravity, jQ; ab- 
solute magnitude, Mboi; radius, R; atmospheric composi- 
tion; etc.). The distances to a number of the closest WDs 
(with D < 100 pc) have been determined unambiguously 
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^ The convention of characterising surface gravity, g, by log g = 
log((;cgs) = log[g/(l cm s~^)] is followed here; the relationship to 
the SI equivalent is log((;si) = log(gcgs) — 2. 



by parallax measure ments (e.g., Bergeron et al.ll200l1 ). and 
the Gaia satellit^l jPerrvman et al.. .2001i) will eventually 
measure pa rallaxes for all Galactic source s brighter than 
r ~ 20 fe.g. jTorres et al.ll2005l : |jordanll2007l '). In most cases, 
however, the characterisation of WDs is necessarily more 
model-dependent, with spectroscopic or photometric data 
used to infer the physical properties of WDs. WD spectra 
encode sufficient information to infer Teg and log g to a rel- 
ative accuracy of < 10 per cent (e.g ., [Bergeron et al .lll992l: 
.McCook fc Si on ,1999: Harris et al'l 12003: Klcinm an et al I 
l2004i : lEisenstein et akluOOSl ), and the a greement between the 
various semi- i ndepe n dent models (e. g., Bergeron et al.lll995l: 
Finlev et all Il997l : iFontaine et al.l I2OOII: iBergeron et all 



2OOII : iBergeron fc Leggettll200 j : iHolberg fc Bergeronll2006D 
gives grounds for confidence in the resultant parameter esti- 
mates. Perhaps the most impressive demonstration of spec- 
troscopic fitting has been the largely auto mated analysis o f 
~ 10* S loan Digital Sky Surve y (SDSS; lYork et all |200GD 
WDs bv lEisenstein et al l (|2006|). This confirmed that most 
WDs have log q ~ 8 (cf. IBergeron et al. 1992: Liebert et al.l 
l2005h . but also that the inferred surface gravities and 
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masses of cool (Teff 5, 10"* K) WDs can be biased high (cf. 



Bergeron et aLlll990l: iKleinman et al.ll2004l : lEisenstein et alj 



200m . Interestin gly, such a bias is not seen in purely pho- 
tometric fits (Engelbrecht fc Koes tcr 2007) . Combined with 
the fact that it is possible to reliably identif y WDs using 
only photometric and astrometric data (e.g., [Harris et al.l 
l2006l : iKilic et al. lEooe), it is of great importance to under- 
stand what constraints can be placed on a WD's properties 
from photometric data alone. 

The most optimistic case is illustrated in Fig. [T] (a), 
in which colour measurements of non- variable point -sources 
generated f r om m ulti-epoch SDSS Stripe 82 data by 
llvezic et al.l (|2007l ) are compared to predicted photom- 
etry generated by G. Fontaine and P. Bergerorlf]. The 
sources shown have relative photometric errors of ~ 0.01 
in all bands, and are mainly WDs, although there is 
some contaminati on at high g—r from quiescent quasars 
ijSesar et al.l 2007i) and unr esolved pairs of white and red 
dwarfs ( Smolcic et al.l |2004 ) . Comparison of the observed 
WD loci with the models in Fig.[T](a) reveals not only an ob- 
vious separation between the H and He WDs, but also clear 
differentiation in Toft . The surface gravity is well constrained 
for H atmosphere WDs with 10* K < Teft < 2 x 10* K (i.e., 
u—g ~ 0.3 and g—r ~ —0.2). Unfortunately, this level of 
characterisation, particularly in log g, relies on Stripe 82's 
near-ideal combination of photometric precision and wave- 
length coverage. Fig.[T](b) shows the results of increasing the 
relative photometric errors of the Stripe 82 data to ~ 0.05: 
the H and He loci are no longer discernable and only the 
broadest constraints could be placed on log g from these 
data. Even more critical than high photometric precision is 
access to near-ultraviolet (NUV) data, as can be seen from 
Fig. [T] (c). In the longer wavelength optical bands WDs' 
colours depend primarily on T^s, with minimal separation 
in log g, or between H and He models. 

However even non-ideal observations of the sort shown 
in Fig. [T] (b) or (c) would allow T^s to be constrained, pro- 
vided only that the WD is uambiguously detected (i.e., with 
relative photometric errors of < 0.1) in at least two pass- 
bands that are blueward of its Rayleigh-Jeans tail. In this 
situation the only other available information comes from 
the flux measurements, which then fix the solid angle, Q,, 
subtended by the source. As 

^ M M 

"«:d2 (1) 

this leads to a well-known degeneracy between distance and 
radius, and hence surface gravity, th at is difficult to break 
with photometric data alone (e.g., iBergeron et al.l 1 19971 . 
I2OOII '). It is, however, possible to use prior knowledge of the 



Galactic WD population - specifically, that most WDs have 
log p ~ 8 - to break this degeneracy. [Harris et al.l (|2006l ) 
adopted this approach in their analysis of a large sample of 
photometric SDSS WDs by assuming log (/ = 8, although 
it is undoubtedly overly-prescriptive to deny any possibil- 
ity that a WD's surface gravity deviates from this fiducial 
value. Nonetheless, this is more realistic than the alterna- 
tive of treating all log g val ues as equally likely a priori; 
this assumption was made bv lVidrih et~aiT(|2007l ). hereafter 
V07, in their analysis of a (differe nt) photometric SDSS WD 
sample. Both [Harris et ahl (I2OO6I ) and V07 identified possi- 
ble halo WDs, although the existence of a strong parameter 
degeneracy means that the identification of outliers on the 
basis of best-fit parameter estimates is likely to yield large 
numbers of spurious candidates. 

The overall implication is that any non-spectroscopic 
WD parameter estimates will depend in part on the fitting 
procedure used, and that prior information about the Galac- 
tic WD population could be as important as the data that 
are specific to the source in question. At least some of these 
issues can be resolved by adopting a more rigorous statis- 
tical approach to estimating WD parameters from photo- 
metric and astrometric data. This rather broad topic will 
be addressed in a series of papers, of which this is the first. 
It deals with the more focussed question: How much can be 
inferred about a WD from photometric data alone? After de- 
scribing the Bayesian formalisms for parameter estimation 
(Section [2| and model selection (Section [Sjl, the WD pop- 
ulation model used as a prior distribution is developed in 
Section [4] After analysing simulated WD observations (Sec- 
tion [5]), these methods are applied to the small samples of 
candidate halo and ultra-cool WDs identified by V07 in Sec- 
tion[6l The conclusions and future directions of this research 
arc summarised in Section [T] 



^ Machine-readable versions of the Fontaine— Bergeron 
H and He atmosphere WD model grids, giving 
broad-band magnitudes for a range of temperatures 
and surface gravities, are currently available from 
http : //www. astro .umontreal . ca/~bergeron/CoolingMod6ls/. 
These are used for all the results presented here as the focus 
is on random statistical uncertainties; hence the the systematic 
variations between Bergeron's mo d els an d the semi-inde p enden t 
models of, e.g . . iBereeron et al. Finlcy ct al. 

[Fontaine et al.l l|2001^ and iHolberg fc Bergeron (,200& ) are not 
explored. 



2 PARAMETER ESTIMATION 

Given photometry of a WD [i.e., measured apparent magni- 
tudes, d — (mi,m2, . . . ,mN^), in each of A^b bands], what 
constraints can be placed on its properties [i.e., the model 
parameters G — {g,Ti,s,D)], assuming its dominant atmo- 
spheric element (H or He) is known? This is a classic param- 
eter estimation problem for which several approaches exist; 
central to all of them is the sampling distribution, or like- 
lihood, which gives the probability of observing the data d 
given a model 0. For the high signal-to-noise ratio measure- 
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Figure 1. Two colour diagrams for a sample of non- variable point— sources in SDSS Stripe 82 from Ivczic et al. (2007). In (a) and (c) the 
sources have relative photometric errors of ~ 0.01 in all bands; in (b) the errors have been inflated to be ~ 0.05. The Fontaine— Bergeron 
model grids for H (red) and He (blue) WDs are also shown, the temperature going from Teg = 110000 K at the bottom-left down to 
Teff ~ 6000 K at the top-right, with the exception of the coolest H models in (c). 



mentfl considered here, the likelihood has the form 
L{d\0) = L{rhi,m2, . . . ,rhNi\g,TcB, D) 

= n (27r)V2A^/ ^P(~x'/2). 



where 



Xb 



2 _ (rhb- [Mt{g,Tcs) + 51og(7:)/10 pc)] 



Arhb 



(2) 



(3) 



Arhj is the estimated noise in the 6th band and Mt(g,Tcfi) 
is given by interpolating between a grid of WD models. The 
effects of the interstellar medium (ISM) are ignored here, 
but could be incorporated into Eq. (|3]) by adding a band- 
dependent extinction term inside the square brackets to ac- 
count for the extra dimming. 

It is at this point that the various methods of param- 
eter estimation diverge. Historically, the most common ap- 
proach has been to use one of several best-fit estimators 
(Section 12. ip . although it is possible to obtain more reliable 
parameter values and uncertainties using Bayesian inference 
(Section [221). 



2.1 Best-fit parameter estimation 

An intuitively appealing and commonly used method of pa- 
rameter estimation is to find the model, 9, that is the best 
fit to the data. This can be found in several different ways, 
including maximizing the likelihood, minimizing x^j or by 
a least squares approach, all of which are equivalent for the 
simple sampling distr ibution given in E q . The lar g e WD 
samples identified bv lKilic et all (j2006l ). lllarris et al.l (|2006l ) 
and V0 7 were analy s ed us ing variants of this approach, al- 
though [HarriT^eF^al] l|2006l ) restricted their investigation to 



In reality the photometric noise includes both a Poisson contri- 
bution due to the finite photon counts and a background uncer- 
tainty that is usually taken to be Gaussian in fiux units. However, 
the WD samples considered here have photometric data accurate 
to a few per cent, sufficient to adopt the approximation that the 
noise is additive and Gaussian in magnitude units. It is further as- 
sumed that there is no covariance between bands: although there 
are measurable inter-band noise correlations in data from a single 
SDSS scan fScranton et al.ii2005iV this effect is far weaker when 
using multi-epoch data. 



models with log g — 8 and V07 exami ned only models on 
Bergeron's grid of pre-computed values. [Harris et af] (|2006l ) 
do not quote uncertainties, and a heuristic procedure to 
obtain uncertainties is used in V07, whereby the best and 
second-best fitting grid points are compared. Given the fun- 
damental reduction in precision attainable when using only 
photometric data to constrain the properties of WDs, it is 
clearly important to investigate this issue more thoroughly. 

2.2 Inferential parameter estimation 

In order to obtain the tightest possible constraints on model 
parameters all the available information must be utilised, 
which implies that any prior knowledge of the model pa- 
rameters must be included (which is also necessary for self- 
consistency; e.g., [jaynes 2003). Applying Bayes's theorem 
to the data, d, yields the posterior probability distribution 



P{0\d) 



n{e)L{d\e) 
E{d) 



(4) 



where n(9) is the prior probability distribution of the pa- 
rameters, L{d\6) is the likelihood, given in Eq. ([2]), and the 
evidence is the integral over the A^p model parameters. 



E{d) = / 7t(0)i(d|0)d6»id6i2 



de. 



(5) 



The evidence plays an important role in model selection 
(Section [3]), but in parameter estimation serves only to nor- 
malise the posterior, and so it is often more convenient to 
work with the unnormalised posterior, defined most simply 



p'{e\d) = n{e)L{d\e) 



(6) 



Formally, the entire posterior distribution is the 
Bayesian answer to any parameter estimation problem, but 
even though it can be written down directly it is not always 
useful in practice. Even with just three model parameters, 
a numerical algorithm is required to determine the location 
and extent of the posterior peak - and the peak is not neces- 
sarily well-defined in the presence of degeneracies. Instead, 
it is far more efficient to explore the paramete r space us- 
ing Markov chai n Monte Carlo (MCMC, e.g., iNeall Il993l ; 
iGilks et al.]|l999l ) techniques. Many different MCMC algo- 
rithms exist, although all have the same aim: to generate a 
sample of model parameter values, {61,82,03, ■■ ■}, drawn 
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from the posterior distribution. For the problem at hand, 
a fairly si mple implementation of the Metropolis- Hastings 
algorithm (jMetropolis et all Il953l : iHastingd Il970l ) is suffi- 
cient, in which the posterior is sampled by a guided random 
walk that tends to explore the high posterior regions whilst 
occasionally exploring the low probability tails of the distri- 
bution. The samples generated from the first stage of this 
process (burn in) must be discarded, as they are determined 
partly by the starting value, but eventually the statistical 
properties of the chain will depend on P' (0\d) alone. Whilst 
there are only heuristic test s to identify these two phases 
of th e MCMC process (e.g., iGelman. A. and Rubin. D. b1 
1 1992! ) a brute force approach can be used to ensure that the 
final chain is well-mixed. 

A Markov chain of samples from the posterior has a 
number of desirable properties, one of the most fundamen- 
tal being that it already samples all of the possible sub- 
sets of model parameters, a task that otherwise would re- 
quire a multi-dimensional marginalisation integral over the 
other parameters. Another advantage of the Monte Carlo 
approach is that the distribution of not only the indepen- 
dent model parameters (i.e., g, T^s and D here) but of any 
derived parameters (e.g., mass, M, or radius, R) can be ob- 
tained simply by calculating the value of interest for each 
sample drawn from the posterior. 

In many situations - in particular that considered here, 
in which chains are generated for each WD in a large sample 
- it is impractical to store the full chains. If the posterior is 
sufficiently localised then it is generally more useful to obtain 
a fiducial parameter estimate, 0, and uncertainty, A0, from 
the chain. As the full Bayesian answer to any parameter es- 
timation problem is the entire posterior distribution, there 
can be no single correct approach; but several algorithms 
give reasonable approximations to single-peaked posteriors. 
The simplest is to take moments of the distribution, thus 
identifying the mean of the distribution as the fiducial es- 
timate 6 and using the variance to obtain AO. Whilst this 
approach is appealingly straightforward, it has the disad- 
vantage of being sensitive to outliers (which can occasion- 
ally be extreme, due to the stochastic nature of the MCMC 
algorithm). A more robust method of generating parame- 
ter estimates and uncertainties is to use the median value 
in each model parameter to construct 6 and then integrate 
the marginalised probability outwards in each direction to 
enclose the desired fraction of the posterior. This approach 
is adopted here, with the quoted uncertainties enclosing the 
central 68 per cent of the marginalised probability in the 
parameter of interest. 



3 MODEL SELECTION 

Section [2] addressed the question of how to estimate the pa- 
rameters of a WD given its dominant atmospheric element, 
but the same data can also be used to determine the likely 
atmospheric composition (or even whether the object is a 
WD at all, although this possibility is not explored here). 
The conventional method of model comparison (Section (ST} 
uses only the quality of the best fit under the different hy- 
potheses, whereas the Bayesian approach marginalises over 
the model parameters and also includes a prior on each 
model (Section 13. 2|) . 



3.1 Best-fit model selection 

The most common approach to model comparison is to use 
the best-fit statistics described in Section r2.1l to characterise 
how well the models fit the data. In the simple case under 
consideration here (i.e., in which the two models have the 
same parameters) the resultant model comparison statistic 
typically reduces to the likelihood ratio. 



exp 



r, (Xmin,H 



Xmiii.HcJ 



(7) 



where the right hand expression assumes the likelihood has 
the form of Eq. In some simple situations this is a suf- 
ficient statistic for model comparison, but the fact that it 
only encodes information about a single point in parameter 
space means that it is of limited utility in general. 

3.2 Inferential model selection 

The central relationship in Bayesian model selection (e.g., 
IJavnedliooal ) is, as in Section [2. 2 1 Bayes's theorem, this time 
applied to the discrete hypotheses that the WD atmosphere 
is predominantly H or He. Given data d and the model pa- 
rameters 6 (which, in this case, are the same for both hy- 
potheses), the probability that the WD atmosphere is H is 



PH(d) 



n-nE-aid) 



nnEn{d) + (1 - nii)EHc{d)' 



(8) 



where tth is the (prior) probability that a WD has a H atmo- 
sphere, the substitution tthc = 1 — tth encodes the assump- 
tion that only H and He atmosphere WDs are considered 
here, and Eu{d) and i?Hc(d) are, as defined in Eq. ([5]), the 
evidence values. 

Whilst the evidence completely characterises how con- 
sistent the data are with each model under consideration, it 
is worth emphasising that the numerical value itself is, like 
the likelihood, not particularly meaningful in isolation. It is 
more useful to scale the evidence relative to its maximum 
possible value (which is a function of the precision of the 
data, as well as the model and priors), but even then a low 
relative value can have two distinct explanations. The most 
intuitive is simply that the data are inconsistent with the 
entire family of models and that that any resultant fit is 
inappropriate; but in most situations there is also the pos- 
sibility that the data are a good fit to a particular subset of 
models which are, a priori, unlikely. 



4 THE WHITE DWARF POPULATION 

In the Bayesian approach to parameter estimation and 
model selection described in Sections [2] and [S] respectively, 
inferences are obtained by combining the available data spe- 
cific to a given WD with information about the Galaxy's 
WD population. If the observations of the source in ques- 
tion are not particularly informative then the prior informa- 
tion on the WD population can be critical, possibly breaking 
the otherwise troublesome parameter degeneracies discussed 
in Section [1] More formally, the WD prior n{g,TcS,D,l,h) 
gives the probability that a WD has a surface gravity of <?, an 
effective temperature of T^s, and a given Galactic position 
(specified by its distance, D and Galactic coordinates, I and 
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6), in the absence of any data specific to the object in ques- 
tion. As such, it is appropriate to use n{g, Tcff , D, I, b) in the 
analysis of any WD, irrespective of how it was discovered 
(e.g., whether it was identified in a flux-limited, volume- 
limited, or more heterogeneous sample). 

In principle, n{g, Tes, D, I, b) should encode a complete 
model of the Galaxy's WD population, giving the relative 
numbers of WDs with different intrinsic properties and po- 
sitions. It should also include the different age and tem- 
perature distributions of WDs drawn from the various sub- 
populations (thin disk; thick disk; halo; etc.). In practice, 
however, Tc{g,Teg, D,l,b) has to be inferred from existing 
surveys and models, and combining the available informa- 
tion to obtain n(g,TcS,D,l,b) is non-trivial. This higher- 
level problem of inferring the properties of the Galaxy's WD 
population is not attempted here, and instead a simpler, 
more heuristic, approach is adopted to obtain a useful ana- 
lytic WD parameter prior. 

The first simplification is to reduce the number of model 
parameters by omitting I and b. This is justified by the 
fact that the angular position of any detected WD is mea- 
sured to a much greater accuracy than the scales on which 
Tc{g,Tes, D,l,b) varies appreciably. In statistical terms the 
astrometric data are considerably more informative than the 
prior; in more astronomical terms, the positions of almost all 
optical sources are known to better than an arcsec, whereas 
no Galactic model would have meaningful structure on such 
small angular scales. However this does not imply that the 
reduced prior, 7i{g, TcS, D), is independent of / and b, a point 
explored further in Section [4.21 

A further plausible simplification is that the reduced 
prior is separable, so that 



n{g,ns,D)=n{g)n{ns)n(D). 



(9) 



The priors on surface gravity and temperature are given in 
Section [TT] the (position-dependent) distance prior is given 
in Section 112] 



4.1 Intrinsic properties 

The most precise and unambiguous constraints on the in- 
trinsic properties of local WDs come from spectroscopic 
studies, with the on- going SDSS campaign descr ibed by 
iKleinman et al.l l|2004) and lEisenstein" et al.l (|2006| ) by far 
the largest such project. The prior distribution of a WD's 
intrinsic properties adopted here is based largely on their 
results. 

Pro bably the most s t raight forward conclusion from the 
study of lEisenstein et al.l ||2006') is that DA WDs (with H at- 
mospheres) outnumber DB WDs (with He atmospheres) by 
~ 10 to 1, with the fractional contribution to the WD pop- 
ulation of other spectral classes being significantly smaller 
again. Thus only H and He models are considered here, and 
the prior probability that a WD has a H atmosphere is taken 
to be Ph = 0.9. 

One of the most important results from the spec- 
troscopic studies of WDs is that they have only a 
small range of surf a ce gravities (e.g., 
iLiebert et al.l l2005l : lEisenstein et al.l 



sented by lEisenstein et al.l l|2006l ) imply that log g for DA 
stars is distributed approximately as 



Tc{g) oc exp 



log 



0.1 



(10) 



Bergeron et"al]|l992l: 



2001 



The fits pre- 



corresponding to a mass range of ~ (0.6 ± 0.1) AIq. There 
appears t o be more scatter for D B than for DA WDs in the 
results of lEisenstein et al] (|2006f ). but it is unclear whether 
this spread is intrinsic. The first major analysis of DBs 
(Beaucha mp et alii 19961 ) implied that they have a narrower 
mass di stribu t ion t han DAs, and the more recent results 
of fVoss et al.1 \200i ) again show a relatively narrow range 
of masses. There is some ambiguity in the surface gravity 
distribut ion due to both the l imitations of the models at 
low TeB (|Eisenstein et al.]|2006 l and references within), and 
the selection effect that all flux-limited samples preferen- 
tially include larger, more luminous, WDs with low log g. 
Indeed, the local vo lume-limited WD sample compiled by 
iHolberg et akl (|2008l ) imphes a log g ~ 8.1, which is in mild 
disagreement with the results obtained by applying a volume 
correc tion to the more rep resentative, but flux-limited, sam- 
ple of iLiebert et al.l l|2005h . It is also possible that there are 
more WDs with \ \ogg- 7.9| > 0.5 than implied by Eq. (fTO)) . 
but the tails of the instrinsic surface gravity distribution are 
not well measured. Given the overall consensus with previ- 
ous studies, Eq. (|10p is adopted here as the informative WD 
surface gravity prior for both DB and DA WDs. 

The present day temperature distribution of WDs is de- 
termined by a combination of their cooling properties and 
the star formation history of the Galaxy, as the lower tem- 
perature white dwarfs have a wide variety of ages. The tem- 
perature distribution of the hottest WDs is simpler because 
they cool very rapidlj^- at a rate of TcB oc T^g'^ accord- 
ing to the models of iMestel (il952) . although the numerical 
cooling curves of the Fontaine-Bergeron models are more 
complex. The implication is that any WD observed with 
TcS > 3 X 10* K formed within the last ~ 10** yr, much 
shorter than the timescale over which the stellar population 
of the Galaxy evolves. Assuming that all WDs are formed 
with the same initial Tes , these simple arguments imply that 
the number density of the hottest, youngest WDs scales as 
dnwD/dToff cx T'J^'^. 

It is non-trivial to estimate the WD temperature dis- 
tr ibution empirical l y, eve n with, e.g., the large WD sample 
of lEisenstein et al.ll|2006h . as the selection effects mentioned 
above in the context of the surface gravity distribution are 
even stronger for temperature due to the L oc T^fg luminos- 
ity dependence. A standard 1/Vmax estimator is only valid 
for a population of objects with uniform spatial distribu- 
tion, whereas most WDs with Tcff > 5 x lO** K would be 
detectable by SDSS at distances of D > 5 x 10^ pc. The tem- 
perature distribution of all but the coolest WDs can only be 
inferred in the context of a full Galactic model (cf. Fig. [2| . 

Fortunately, optical and near-UV photometry is gener- 
ally sufficient to unambiguously constrain a WD's effective 
temperature, at least up to Toff ~ 6 x 10* K, so the pa- 
rameter estimation results described in Section[2]should not 
be greatly affected by this uncertainty in the prior. Above 
this temperature, these bands cover only the Rayleigh- Jeans 
tail of a WD's spectrum, and the measured colours can 
only provide a lower limit on the effective temperature of 
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Teft <; 6 X 10'' K. In a second fortunate coincidence it is in 
just this regime that the temperature distribution inferred 
from the above cooling argument is likely be reasonably ac- 
curate, the lower limit on Toff implied by the data being 
complemented by the sharp T~^^^ dependence of the prior. 
One of the eventual aims of this research is to apply the 
statistical approach described in Section [2]2] to constraining 
the properties of the WD population, but for the moment 
the temperature prior is taken to be 



7r(refr) (X 



-7/2 



(11) 



4.2 Galactic distribution 



As WDs are the evolutionary endpoints of MS stars it is 
expected that they have a similar distribution in Galactic 
position and velocity, mostly being located in - and rotating 
around the Galactic Centre with - the stellar disk(s). This 
is demonstrably the case for most known WDs, although 
the fact that WDs are so long-lived leads to the intriguing 
possibility of a significant relic population of cool, old, halo 
WDs fe.g..lChabrieij|l999l : lAlcock et al]|2000l : iHodgkin et al] 



l2000l : llbata et al.ll2000l V There have been claims that such a 
population has been detected ([Opp onheim or ct al. 2001a b), 
although several studies have disputed the inferred sp ace 
density of this population (jSalim et al.ll2004l : lReidll2005l and 
references therein). 

The spatial distribution of WDs adopted here includes 
both thin and thick stellar disks, along w ith the halo , and i s 
based on the Galactic model described bv lJuric et al.l (|2008l ) . 
In cylindrical Galactic coordinates R, Z and $, the number 
density of stars per unit volume is taken to scale as 



nMs{R,Z, $) 



exp 



R 



+./'tiiick exp ( — 



-I- ./halo 



exp 



R 



exp 



ithick / \ ^^thick 

[i?' + (^/ghalo)'+flLc]''" 



,(12) 



where /thick = 0.06, /halo = 6 x 10 ^, ithin = i/thick = 
3500 pc, Lhaio = 8500 pc, i?corc = 1000 pc, i/thin = 260 pc, 
-ffthick = 1000 pc, (jhaio = 0.64 and 77haio = 2. This den- 
sity can then be expressed in terms of heliocentric Galactic 
coordinates, defined implicitly by 

R{D,l,b) = [D^cos^(6) - 2RqD cos{b) cos{l) + R%Y^^ , 

Z{D,l,b) = Dsm{b), (13) 

I, b) = arctan [D cos(6) cos(i) — Rq,D cos(&) sin(/)] , 

where Rq — 8000 pc is the fiducial distance from the Sun 
to the Galactic centre. The fact that the Sun is actu- 
ally 'above' the Galac tic plane with Zq = (24.2 ± 1.7) pc 
(|Mafz-Apellanij l200lh is unimportant here as the stellar 
density varies appreciably only over much greater scales than 
this. The distribution in D, I and b that results from substi- 
tuting these definitions into Eq. (|12|l is rather cumbersome, 
and so not given explicitly here. As argued above, the prior 
distribution in I and b is unimportant, whereas implied dis- 
tance distribution has the potential to affect the parameter 



estimates for a WD, especially given the solid angle degen- 
eracy discussed in Section [1] The distance prior is obtained 
by multiplying Eq. (|12p by the volume element (ex D^) to 
give 

Tv{D\l, b) oc hms IR{D, I, b), Z{D, I, b), $(D, I, b)] , (14) 

which is shown in in Fig. [2] for several values of I and b. 

The low-D behaviour is determined by the volume ele- 
ment; beyond ~ 100 pc the Galactic direction is more impor- 
tant, with the prior rising steeply towards the Galactic cen- 
tre, but dropping off sharply when looking out of the plane. 
In all cases the high-D tail is dominated by the halo com- 
ponent, which is possibly an underestimate for WDs, but 
which should also fall off faster than this power-law form 
given m Eq. dHJ for _D > 10^ pc. A more realistic model 
would include a sharper fall off in the halo density, but this 
is irrelevant for the relatively shallow surveys (which can 
detect WDs only out to D ~ 10* pc) considered here. 

One immediate, if somewhat counter-intuitive, impli- 
cation of the directional dependence of the distance prior 
is that the parameters inferred for two WDs with identi- 
cal photometry could be quite different if they are on well 
separated lines-of-sight. Consider, for example, a very faint 
WD with measured colours that, if it is assumed to have 
log g — 8, imply a distance of D ~ 10'* pc. Observed towards 
the Galactic centre this would be completely consistent with 
prior expectation (cf. the I = deg; 6 = deg curve in Fig. [2] 
a). Observed out of the Galactic plane this would be very 
surprising, as there are expected to be very few WDs at this 
distance (which is much greater than the disk scale heights) 
above or below the Galactic plane, and it would be much 
more likely that the object is a closer, and correspondingly 
more massive, WD with log g — 9. Most WDs identified in 
the current generation of wide-field surveys are closer than 
this (i.e., local) and so the variation in the Galactic den- 
sity field is less relevant, but it still has some influence in 
identifying potential halo WDs. 



4.3 Discussion 

Whilst the prior given by combining Eqs (I10|l , and (|14|l 
is a reasonable approximation to the Galaxy's WD popula- 
tion, it is incomplete in several ways. The separation of the 
intrinsic distributions from the positional distribution prob- 
ably leads to the numbers of thick disk and halo WDs being 
underestimated (as these stellar populations are older than 
that of the thin disk). This model is also equivalent to as- 
suming that the different Galactic WD populations have the 
same age and temperature distributions, which is probably 
unrealistic, ev en though the re is considerable uncertainty in 
this area (e.g. . iHansenl 200ll ) . However these approximations 
are only important here to the degree that they affect the 
eventual WD parameter estimates. For photometric data of 
the quality considered in Section [5] and [6] the most impor- 
tant role of the prior is to encode the flducial surface gravity 
value and the relative numbers of thin disk, thick disk and 
halo WDs, for which Eqs dTO]) and HH) suffice. 
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Figure 2. The prior distributions of D that are used here, based on the Galactic model of Juric et al. (2008), and described in Section|4] 
The prior along four fiducial lines-of-sight is shown in (a). The prior along the lines— of— sight to the middle and both ends of SDSS 
Stripe 82 (from which the V07 sample is drawn) is shown in (b). The different lines— of— sight are labelled: by Galactic longitude, I and 
latitude, b, in (a); and by right ascension (as Stripe 82 is at declination 5 = 0) in (b). Note that the vertical axes differ, and that the 
distributions are not normalised, instead being matched to have the same low-D dependence. 
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Figure 3. Marginalised posterior distributions of the model parameters for a simulated WD with T^g = 10200 K and logg = 7.9 at 
D = 179 pc, observed in the g, r, i, z, J, H, K bands with photometric noise of Am = 0.01. For both the non-informative prior (grey) 
and the informative prior (orange) the coloured regions enclose 68 per cent and 95 per cent of the posterior probability in these two 
parameters. Left: the constraints on the independent model parameters T^g and logg are shown, along with the points on the Fontaine- 
Bergeron model grid (crosses). Right: the degeneracy for black body models (black line) is highlighted by showing constraints on the 
WD's distance, D, and radius, characterised by log g — log(Af/MQ). In both panels the true model is indicated by the black spot. 



5 SIMULATIONS 



The implication of the parameter inference and model selec- 
tion methods (described in Sections[2]and[3} for WD surveys 
can be further understood by applying them to simulated 
data. A number of representative cases are examined here, 
although care must be taken in interpreting the Bayesian 
results obtained using the Galactic WD prior given in Sec- 
tion|4l as the simulated WDs have not been drawn from this 
distribution. 



solute AB magnitudes in fou][f| SDSS bands {g, r, i and z) 
and three NIR bands (J, H and K), which were converted 
to apparent magnitudes by placing the WD at distance D. 
Finally, photometric noise (of the same level. Am = 0.01, in 
each band, as appropriate to the SDSS Stripe 82 data) was 
added to give the simulated data g, f, i, z, J, H , and K. 
The data were then analysed to give parameter constraints 
(Section 15. ip and also to test the degree to which H and He 
WDs could be distinguished (Section |5.2|) . 



In each simulation, a surface gravity, g, and an effec- 
tive temperature, T^g, were chosen from within the grid of 
H atmosphere models provided by Bergeron. This gave ab- 



^ The u band was not included in the simulations in order to best 
replicate the analysis of V07 described in Section [6] 
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5.1 Parameter estimation 

WD parameter estimates were obtained from the simula- 
tions in several different ways. Both the best-fit model and 
the best-fit grid value were calculated as described in Sec- 
tion 12. f I and the full posterior distributions of the model 
parameters were found using the MCMC algorithm outlined 
in Section [2.21 In the latter case the analysis was performed 
both with a non-informative prior (uniform in the model pa- 
rameters) and with the full Galactic WD population prior 
given in Eq. ([9j . 

The first case considered is that of a fiducial WD with 
log g = 7.9 and TcS = 10200 K, results for which are shown 
in Fig. 131 The most striking feature of this simulation is 
the strong degeneracy apparent in Fig. [3] (b), which fol- 
lows the form given in Eq. ([T]). The temperature is well 
constrained for both the non-informative and informative 
priors, although the weak correlation between T^s and log g 
means the constraints are better in the latter case. 

This example also illustrates an important point about 
the identification of outliers, in particular fast-moving halo 
WDs. If this WD was analysed using the grid-based best- 
fit method described in Section [2.11 it would be inferred to 
have logg = 7.5 and Tcs — lOQQO K, simply because this 
model grid point is closest to the locus of high likelihood. 
The resultant distance estimate would be biased, high in 
this case, by a factor of D/D ~ 1.2. The tangential velocity 
would also be over-estimated, with the possible result that 
the WD would be identified as a halo object if it had a high 
proper motion. A more extreme form of this effect explains 
several of the candidate halo WDs analysed in Section [6.11 
It is also an illustration of the general principle that it is 
no more valid to infer that a WD is dynamically interest- 
ing due a large best-fit distance, than it is to infer that the 
same star is intrinsically unusual by restricting its distance 
through the application of a prior on its motion. The two 
possibilities should be balanced by giving the appropriate 
(not necessarily equal) weight to all the available informa- 
tion. 

The second example, shown in Fig. |4l is of a T^h — 
10200 K WD with an unusually high surface gravity of 
log g = 8.6. The results are similar to those for the first ex- 
ample, except for the fact that the true model is well away 
from the peak of the posterior obtained with the informative 
prior. This comes about as the prior based on the Galactic 
distribution of WDs (with log g ~ 8) is inappropriate to 
the case of an arbitrarily chosen example. If a real WD was 
observed to have these colours and the log g prior given in 
Eq. (|10p is correct, this would be the appropriate inference. 
In that case the small number of WDs with true colours 
that match the data would be out-numbered by the fiducial 
WDs with log g — 8 that were randomly scattered to these 
observed values. 

The resuhs for a hot WD with T^a = 61000 K and 
log (? =8.1 are shown in Fig. [S] The principal difference 
here is that the WD's emission peaks at a wavelength of 
A ~ 0.05 ^.m, significantly blueward of the SDSS g band, 
which has a peak response at A ~ 0.5 \im. As a result it 
is only the source's Rayleigh-Jeans tail which is measured, 
and any temperature of > 6 x 10* K is consistent with the 
data. Both the best-fit analysis and the Bayesian calculation 
with the non-informative priors give essentially arbitrary pa- 
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Figure 4. Same as Fig. \3\ but for a simulated WD with log g - 
8.6 and T^a = 10200 K at D = 179 pc. 
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Figure 5. Same as Fig. |3] but for a simulated WD with T^ft 
61000 K and log g = 8.1 at D = 463 pc. 



rameter estimates, although the latter does at legist encode 
the fact that the uncertainties are large. However when the 
WD population model from Section |3] is applied, the lower 
bound on Tcff implied by the data combines with the prior 
observation that high- Toft WDs are rare to give a far more 
precise estimate. This is a particularly good illustration of 
the equal role that the prior and likelihood have in Bayesian 
parameter estimation: in isolation they only provide upper 
and lower bounds, respectively; together they combine the 
available information to provide the best possible estimate 
of the WD's temperature. 

Results for the final example, a cool WD with Tcff = 
4150 K and log g = 7.8, are shown in Fig. [B] At these tem- 
peratures the WD colours have a much greater sensitivity to 
g, and so the data are almost completely informative, with 
the result that the parameter estimates are insensitive to the 
nature of the prior, or even to the choice of analysis method. 
The solid angle degeneracy is still apparent, although the 
sparsity of the grid (relative to the parameter uncertainties) 
is such that the best fit from a non-interpolative grid-based 
search is likely to be outside the formal confidence region in 
this temperature range. 
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Figure 6. Same as Fig. |3] but for a simulated WD with T^f{ 
4150 K and log g = 7.8 at D = 84 pc. 



5.2 Model selection 

The model selection techniques described in Section [3. 2 1 can 
be used in a variety of ways; only one illustrative example is 
considered here, which is to assess the importance of u band 
photometry to WD parameter estimation. 

As di scussed in Se ction [1] the SDSS Stripe 82 data pre- 
sented bv llvezic et al.l (2007) suggests that per cent level u, 
g and r band photometry should be sufficient to both de- 
termine the surface gravity of H atmosphere WDs to rea- 
sonable precision, and to distinguish between WDs with 
H and He atmospheres, at least in the temperature range 
10000 K < Tcff < 18000 K. Both these possibilities were ex- 
plored quantitatively by simulating observations of a WD 
with log g = 8.0 and Tcb = 13000 K, using the methodology 
described above, but using all the SDSS bands (i.e., u, g, r, 
i and z). 

The results for parameter estimation, obtained with the 
non- informative priors described in Section |4j are shown in 
Fig. [7] This confirms that accurate optical and near-UV pho- 
tometry is sufficient to estimate log 5 to within ~ 0.1, al- 
though a second simulation in only g, r, i and z also makes 
it clear that the u band is critical to doing so. The loss of the 
u band also results in increased uncertainties in T^b as the 
remaining bands do not probe the Wien part of the WD's 
spectrum as effectively. These points are particularly rele- 
vant to the sample of candidate halo WDs analysed by V07 
and reanalysed here in Section |6l 

Given that the separation between the H and He models 
in Fig. [l]is far greater than that between the various H mod- 
els with Tcft — 13000 K, it is clear that photometric data 
should allow the two classes of models to be distinguished. 
The evidence, defined in Eq. ([5}, was calculated for H and 
He models using the informative priors defined in Eq. (|9]), 
both with and without u band data. With u measurements 
the results are as decisive as expected: En ~ 0.052 and 
Eue — 0.0, implying almost no chance that this WD has a 
He atmosphere, with Pho — 10~^. However without the u 
band data the situation is not quite as clear-cut: Eh — 7.8 
and Ehc — 1-4, implying Phc — 0.02. 




T.,,/10^ K 

Figure 7. Marginalised posterior distributions of the model pa- 
rameters for a simulated WD with T^ff = 13000 K and log g = 8.0 
at D = 200 pc, observed with photometric noise of Am = 0.01 
in u, g, r, i and z (shaded contours) and g, r, i and 2 only (open 
contours) The contours enclose 68 per cent and 95 per cent of the 
posterior probability (obtained using non-informative priors) on 
the independent model parameters T^g and log g. The true model 
is indicated by the black spot. 



6 OUTLIERS 

One of the main utilities of the statistical methods described 
in Sections[2]and[3]is the assessment and characterisation of 
putative outliers from the population under consideration. 
An example of a problem for which this approach is par- 
ticularly well suited is that of identifyi ng halo WDs. This 
subjec t has provoke d much debate fe.g.. lOppenheimer et al.l 
l2001al lbl: lReidll20()5l ) , in part because of the variety of statis- 
tical methods used. 

One of the largest systematically selected samples of 
candidate halo objects is that generated by V07, who 
also identified a number of candidate ultra-cool, and 
hence potentially old, WDs. These samples were extracted 
from the SDSS Lig ht-Motion Curve Catalogue (LMCC; 
iBramich et al.l [iooil ). which provides data for ~ 4 x 10® 
sources in SDSS Stripe 82 that were observed at ~ 20 dif- 
ferent epochs spanning ~ 7 yr. As such, the LMCC has 
improved astrometry and photometry relative to the main 
SDSS survey, with uncertainties of Aft ~ 5 milliarcsec yr~^ 
and Af ~ 0.01 at f ~ 20. The photometry in the g, i and z 
bands is of similar quality, although there were fewer good u 
band observations, necessitating the exclusion of the u band 
from the V07 analysis, a significant loss in the context of 
WD studies (cf. Section [5.2p . V07 were able to produce an 
exquisite reduced proper motion (RPM) diagram in which 
the WD population was very clearly separated from the MS 
stars, and the resultant sample of 1049 WDs is probably 
contaminated at no more than the ~ 2 per cent level (cf. 
iKihc et al.|[2006l : [Harris et al.ll2006l ). 

V07 then cross-matched the WD sample with 
the UKIRT I nfrare d Deep Sky Survey (UKIDSS; 
[Lawrence et al.l [2003), which yielded Y, J, H and K 
band photometry for about half the sources (as the 
UKIDSS observations in Stripe 82 were only partially 
complete at the time). The addition of NIR photometry 
is particularly valuable in the analysis of the coolest WDs 
due to the observed depleted fiux at these wavelengths 
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l|Hodgkin et aLll2000l ). although the Y band data were not 
used by V07 as the Fontaine-Bergeron model grids do not 
include this filter. 

The next stage in V07's analysis was to estimate the 
WDs' parameters from the stacked g, r, i and z band 
SDSS photometry and, where available, the J, H and K 
UKIDSS measurements. This was done by finding the best 
fit from a discrete set of models defined by combining the 
Fontaine-Bergeron grid with a uniform sampling (multi- 
ples of 5 pc) in distance. This yielded parameter estimates 
logg, Teff and D which, when combined with the observed 
proper motion, fi, yielded a tangential velocity estimate of 
Vtan = Dp.. Uncertainties were characterised by finding the 
location of the second-best fitting grid point, although these 
values were not explicitly reported. V07 then identified 34 
candidate halo WDs with estimated tangential velocities of 
Vtan ^ 160 km s~^, as well as 24 candidate ultra-cool WDs 
with best-fit temperatures of T^s ^ 4000 K. V07 took care 
to point out some of the systematic effects, including the 
solid angle degeneracy, that might affect the candidate se- 
lection, and one of the motivations for this reanalysis is to 
quantify some of these biases. 

6.1 Candidate halo white dwarfs 

The 34 candidate halo WDs identified by V07 were reanal- 
ysed using the Bayesian techniques described in Sections [2] 
and [31 giving parameter estimates and evidence values for 
both H and He atmosphere models, and for both the non- 
informative and informative priors defined in Section 3] The 
results are presented in Tables[T]and[51 The Bayesian param- 
eter fits are compared to those obtained by V07 in Figs. |8l 
[Qland fTOl The model comparison results, used to distinguish 
between the H and He models, are presented in Fig. 1111 

Fig. [8] shows the expected result that the Bayesian ef- 
fective temperature estimates match those of V07. The es- 
timates are insensitive to the prior, as the well-measured 
SDSS (and UKIDSS) colours unambiguously constrain Tcs. 
If the data are sufficiently informative then the details of 
the parameter estimation method are secondary, and that 
is the situation here. The clear agreement between the two 
estimates notwithstanding, the Bayesian values are system- 
atically higher; this is an artefact of the much stronger dis- 
tance bias discussed below, which has a small effect here 
because there is a slight anti-correlation between T^s and D 
for WDs with Tcff ~ 10'' K. 

The Bayesian distance estimates for the halo candidates 
differ considerably from those given by V07, as can be seen 
from Fig. [9] For both the non-informative and informative 
priors the Bayesian estimates are lower, for what turn out 
to be related (but distinct) reasons. 

In the case of non-informative priors it is because the 
V07 parameter search is performed on a discrete grid of 
models. The grid point with the lowest is unlikely to be 
the true best-fit, even if it might be close. In many problems 
this would merely result in a small error in the parameter es- 
timates on a scale comparable to the grid spacing. However, 
the strong degeneracy in the WD models exacerbates this 
problem, as the locus of models with good fits can 'thread' 
the model grid, as shown in Figs. [31 [H and [SI This explains 
why the grid-based estimates can differ significantly from 
the true best fit distance, but not why they are systemati- 
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Figure 8. The estimated effective temperature, Tefli of the 34 
candidate halo WDs as found by V07, and using the Bayesian 
methods described here (for the H atmosphere models with the 
Galactic model priors). 



cally higher. The bias seen in Fig. [9] can be traced back to 
the fact that the 34 candidate halo WDs were not selected 
randomly, but by virtue of having large Vtan and hence high 
D. Thus WDs for which the error due to the grid-based fit- 
ting randomly leads to an over-estimate of D will preferen- 
tially be selected into a sample of objects with high inferred 
tangential velocities, explaining the observed bias. 

The difference between the Bayesian distance estimates 
obtained using the informative prior and those of V07 are 
even greater, as shown in Fig. [9](b). The selection bias de- 
scribed above is still present, but the strongly peaked \ogg 
prior given in Eq. (|10p is the dominant factor. With all the 
models along the degeneracy being comparably good fits to 
the data, the informative prior results in the posterior being 
strongly peaked at, in this case, log g ~ 7.9. This effectively 
breaks the degeneracy, and the tight limits on log g from the 
prior then translate into a comparably tight constraint on 
D. The obvious locus of points with D ~ 0.6 -Dgrid are those 
for which the best fit grid point has log^ = 7.0 and the 
Bayesian models have the prior-selected value of log g ~ 7.9; 
the distance ratio implied by Eq. ((1} is 0.59. There are also 
a few points with D ~ 0.85Dgrid; these are the WDs for 
which the best-fit grid point has \ogg = 7.5. Whilst it is 
not necessarily unexpected that a warm, and hence young, 
halo WD would have low mass and log g, such objects are 
expected to extremely rare. A more probable explanation is 
that these are members of the much greater thin disk and 
thick disk populations for which the grid-based fitting pro- 
cedures happened to yield the the twin abnormalities of high 
D and low log g. (For a more detailed discussion of this idea 
see Section [6.1. in The tightness of the loci in Fig. [9] is also 
rather striking, and suggests that the informative log g prior 
defined in Eq. (|10p is overly prescriptive; the true log g dis- 
tribution almost certainly has broader, non-Gaussian wings. 
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Table 1. Estimated H atmosphere model parameters of the 34 V07 candidate halo WDs, obtained using the 
informative prior defined in Eq. ifo) . 
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88 


± 





10 





55 


+ 





05 


(3 


09 


± 





30' 


X 


10* 


496 ± 27 


115 ± 13 


SDSS 


J005906 


78+001725 


2 


(1 


07 


± 





01^ 


X 


10* 


7 


88 


± 





10 





54 


+ 





05 


(4 


48 


± 


Q 


42' 


X 


10* 


294 ± 16 


143 ± 10 


SDSS 


J010129 


81-003041 


7 


(1 


28 


± 





03^ 


X 


10* 


7 


88 


± 





10 





55 


+ 





05 


(2 


76 


± 


Q 


28' 


X 


10* 


565 ± 31 


107 ± 15 


SDSS 


J010207 


24-003259 


7 


(1 


08 


± 





01^ 


X 


10* 


7 


89 


± 





10 





55 


+ 





05 


(4 


40 


± 





41) 


X 


10* 


191 ± 11 


334 ± 19 


SDSS 


J010225 


13-005458 


4 


(1 


43 


± 





03^ 


X 


10* 


7 


88 


± 





10 





56 


+ 





05 


(1 


97 


± 





22) 


X 


10* 


429 ± 24 


216 ± 15 


SDSS 


J014247 


10+005228 


4 


(1 


27 


± 





02* 


X 


10* 


7 


87 


± 





10 





55 


+ 





05 


(2 


81 


± 





28) 


X 


10* 


408 ± 22 


114 ± 10 


SDSS 


J015227 


57-002421 


1 


(9 


21 


± 





11) 


X 


10^ 


7 


89 


± 





10 





54 


± 





05 


(6 


69 


± 





65) 


X 


10* 


351 ± 21 


118 ± 11 


SDSS 


J020241 


81-005743 





(7 


59 


± 





07) 


X 


10^ 


7 


89 


± 





10 





54 


± 





05 


(1 


10 


± 





11) 


X 


lO** 


226 ± 13 


113 ± 8 


SDSS 


J020729 


85+000637 


6 


(5 


64 


± 





03) 


X 


10^ 


7 


88 


± 





10 





53 


± 





05 


(2 


50 


± 





31) 


X 


109 


163 ± 9 


115 ± 8 


SDSS 


J024529 


69-004229 


8 


(1 


42 


± 





03) 


X 


10* 


7 


87 


± 





10 





55 


± 





05 


(1 


98 


± 





22) 


X 


10* 


568 ± 32 


137 ± 13 


SDSS 


J024837 


53-003123 


9 


(8 


80 


± 





07) 


X 


10^ 


7 


89 


± 





10 





54 


± 





05 


(7 


54 


± 





73) 


X 


10* 


216 ± 13 


328 ± 20 


SDSS 


J025325 


83-002751 


5 


(1 


43 


± 





03) 


X 


10* 


7 


80 


± 





10 





52 


± 





05 


(1 


77 


± 





20) 


X 


10* 


265 ± 15 


125 ± 8 


SDSS 


J025531 


00-005552 


8 


(1 


13 


± 





02) 


X 


10* 


7 


88 


± 





10 





55 


± 





05 


(3 


88 


± 





38) 


X 


10* 


448 ± 25 


149 ± 12 


SDSS 


J030433 


61-002733 


2 


(6 


46 


± 





05) 


X 


10^ 


7 


89 


± 





10 





54 


± 





05 


(1 


65 


± 





18) 


X 


109 


283 ± 16 


130 ± 12 


SDSS 


J211928 


44-002632 


9 


(1 


52 


± 





04) 


X 


10* 


7 


87 


± 





10 





55 


± 





05 


(1 


60 


± 





20) 


X 


10* 


461 ± 27 


103 ± 12 


SDSS 


J213641 


39+010504 


9 


(1 


16 


± 





02) 


X 


10* 


7 


88 


± 





10 





55 


± 





05 


(3 


55 


± 





34) 


X 


10* 


454 ± 25 


97 ± 14 


SDSS 


J215138 


09+003222 


3 


(7 


22 


± 





06) 


X 


10^ 


7 


89 


± 





10 





54 


± 





05 


(1 


26 


± 





13) 


X 


109 


247 ± 15 


124 ± 10 


SDSS 


J223808 


18+003247 


9 


(6 


53 


± 





05) 


X 


10^ 


7 


88 


± 





10 





53 


± 





05 


(1 


59 


± 





18) 


X 


109 


213 ± 12 


289 ± 18 


SDSS 


J223815 


97-011336 


9 


(7 


61 


± 





07) 


X 


10^ 


7 


88 


± 





10 





53 


± 





05 


(1 


08 


± 





11) 


X 


109 


328 ± 20 


110 ± 13 


SDSS 


J230534 


79-010225 


2 


(1 


04 


± 





01) 


X 


10* 


7 


89 


± 





10 





55 


± 





05 


(4 


88 


± 





46) 


X 


10* 


431 ± 24 


145 ± 13 


SDSS 


J231626 


98+004607 





(4 


88 


± 





03) 


X 


10^ 


7 


89 


± 





10 





52 


± 





05 


(4 


99 


± 





94) 


X 


109 


130 ± 7 


100 ± 6 


SDSS 


J233227 


63-010713 


8 


(1 


18 


± 





02) 


X 


10* 


7 


88 


± 





10 





55 


± 





05 


(3 


47 


± 





34) 


X 


10* 


379 ± 21 


97 ± 9 


SDSS 


J233817 


06-005720 


1 


(9 


75 


± 





12) 


X 


10^ 


7 


88 


± 





10 





54 


± 





05 


(5 


75 


± 





55) 


X 


10* 


440 ± 26 


177 ± 16 


SDSS 


J234110 


13+003259 


5 


(1 


21 


± 





02) 


X 


10* 


7 


87 


± 





10 





54 


± 





05 


(3 


18 


± 





32) 


X 


10* 


342 ± 19 


177 ± 12 


SDSS 


J235138 


85+002716 


9 


(7 


08 


± 





06) 


X 


10^ 


7 


89 


± 





10 





54 


± 





05 


(1 


32 


± 





14) 


X 


109 


321 ± 19 


146 ± 15 
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Figure 9. The estimated distance, _D, to the 34 candidate halo WDs as calculated by V07 and obtained using the Bayesian method 
described here [for H atmosphere models, with the non-informative priors in (a) and the Galactic model priors in (b)]. The two dotted 
lines show the loci D = 0.6-Dgrid and D = 0.85-Dgrid, expected for WDs with unconstrained surface gravities for which V07 estimate 
logg = 7.0 and logg = 7.5, respectively. 
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Table 2. Estimated He atmosphere model parameters of the 34 V07 candidate halo WDs, obtained using the 
informative prior defined in Eq. l[9]l. 



object T„ff logs M T D Vtan 

(K) (Mg) (yr) (pc) (kms-1) 



SDSS 


J000244 


03+010945 


8 


(1 


12 







02) 


X 


10*" 


7 


89 







10 





53 


± 





05 


(4 


29 


± 





52) 


X 


10*^ 


436 




27 


105 ± 12 


SDSS 


J000557 


24+001833 


2 


(8 


87 







08) 


X 


10"^ 


7 


89 







10 





53 







05 


(8 


02 







94) 


X 


10^ 


172 




10 


171 ± 11 


SDSS 


J001306 


23+005506 


4 


(1 


07 







01) 


X 


10* 


7 


89 







10 





53 







05 


(4 


92 







59) 


X 


10^ 


306 




19 


152 ± 11 


SDSS 


J001518 


33_|_010549 


2 


n 


27 


+ 





02^ 




lO'' 


7 


89 


+ 





10 





54 


+ 





05 


(3 


09 


+ 





40) 




10® 


472 




30 


154 ± 15 


SDSS 


J001838 


54-1-005943 


5 


(1 


96 







07) 


X 


lO"* 


7 


90 







10 





55 







04 


(7 


42 




1 


31) 


X 


107 


678 


± 


44 


109 ± 18 


SDSS 


J002951 


64+005623 


6 


(1 


19 


± 





02) 


X 


10* 


7 


89 


± 





10 





54 


± 





05 


(3 


69 


± 





47) 


X 


10'* 


535 


± 


34 


221 ± 21 


SDSS 


J003054 


06+001115 


6 


(7 


84 







08) 


X 


10"^ 


7 


90 







10 





53 







05 


(1 


11 







14) 


X 


lO** 


349 




21 


159 ± 16 


SDSS 


J003730 


58—001657 


8 


(8 


75 







10) 


X 


10"^ 


7 


89 







10 





53 







05 


(8 


34 


± 





99; 


X 


10*^ 


297 




18 


106 ± 9 


SDSS 


J003813 


53-000128 


8 


(1 


25 


± 





02) 


X 


10* 


7 


89 


± 





10 





53 


± 





05 


(3 


18 


± 





4o; 


X 


10* 


308 


± 


19 


219 ± 15 


SDSS 


J004214 


88-1-001135 


7 


(1 


24 







02) 


X 


10* 


7 


89 







10 





53 







05 


(3 


25 


± 





43) 


X 


10® 


486 




31 


113 ± 14 


SDSS 


J005906 


78+001725 


2 


(1 


08 







01) 


X 


10* 


7 


89 







10 





53 







05 


(4 


73 


± 





57) 


X 


10® 


284 


± 


18 


138 ± 10 


SDSS 


J010129 


81 — 003041 


7 


(1 


29 







02) 


X 


10* 


7 


89 







10 





54 


± 





05 


(2 


91 







39) 


X 


10® 


559 


± 


35 


106 ± 15 


SDSS 


J010207 


24—003259 


7 


(1 


10 


± 





01) 


X 


10* 


7 


89 


± 





10 





53 


± 





05 


(4 


54 


± 





54) 


X 


10® 


186 


± 


11 


327 ± 20 


SDSS 


J010225 


13—005458 


4 


(1 


47 


± 





03) 


X 


10* 


7 


89 


± 





10 





54 


± 





05 


(1 


97 


± 





29) 


X 


10® 


442 


± 


29 


222 ± 17 


SDSS 


J014247 


10+005228 


4 


(1 


29 


± 





02) 


X 


10* 


7 


89 


± 





10 





54 


± 





05 


(2 


91 


± 





38) 


X 


10® 


405 


± 


25 


113 ± 10 


SDSS 


J015227 


57-002421 


1 


(9 


13 


± 





11) 


X 


10^ 


7 


89 


± 





10 





53 


± 





05 


(7 


46 


± 





88) 


X 


10® 


340 


± 


21 


115 ± 11 


SDSS 


J020241 


81-005743 





(7 


46 


± 





06) 


X 


10^ 


7 


89 


± 





10 





52 


± 





05 


(1 


25 


± 





15) 


X 


109 


222 


± 


13 


112 ± 8 


SDSS 


J020729 


85+000637 


6 


(5 


67 


± 





03) 


X 


10^ 


7 


93 


± 





10 





54 


± 





05 


(3 


06 


± 





46) 


X 


109 


164 


± 


10 


116 ± 8 


SDSS 


J024529 


69-004229 


8 


(1 


45 


± 





03) 


X 


10* 


7 


90 


± 





10 





54 


± 





05 


(2 


08 


± 





30) 


X 


10® 


573 


± 


37 


139 ± 14 


SDSS 


J024837 


53-003123 


9 


(8 


67 


± 





08) 


X 


10^ 


7 


90 


± 





09 





53 


± 





05 


(8 


67 


± 


1 


01) 


X 


10® 


208 


± 


12 


316 ± 19 


SDSS 


J025325 


83-002751 


5 


(1 


49 


± 





03) 


X 


10* 


7 


89 


± 





10 





54 


± 





05 


(1 


88 


± 





27) 


X 


10® 


263 


± 


17 


125 ± 9 


SDSS 


J025531 


00-005552 


8 


(1 


14 


± 





02) 


X 


10* 


7 


89 


± 





10 





53 


± 





05 


(4 


09 


± 





50) 


X 


10® 


432 


± 


27 


143 ± 12 


SDSS 


J030433 


61-002733 


2 


(6 


40 


± 





04) 


X 


10^ 


7 


91 


± 





10 





53 


± 





05 


(1 


88 


± 





24) 


X 


109 


280 


± 


17 


129 ± 12 


SDSS 


J211928 


44-002632 


9 


(1 


56 


± 





04) 


X 


10* 


7 


90 


± 





10 





55 


± 





05 


(1 


64 


± 





25) 


X 


10® 


474 


± 


32 


106 ± 13 


SDSS 


J213641 


39+010504 


9 


(1 


18 


± 





02) 


X 


10* 


7 


89 


± 





10 





53 


± 





05 


(3 


75 


± 





47) 


X 


10® 


438 


± 


28 


94 ± 14 


SDSS 


J215138 


09+003222 


3 


(7 


12 


± 





06) 


X 


10^ 


7 


90 


± 





10 





53 


± 





05 


(1 


41 


± 





17) 


X 


109 


244 


± 


15 


123 ± 10 


SDSS 


J223808 


18+003247 


9 


(6 


43 


± 





04) 


X 


10^ 


7 


93 


± 





10 





54 


± 





05 


(1 


89 


± 





24) 


X 


109 


206 


± 


12 


279 ± 18 


SDSS 


J223815 


97-011336 


9 


(7 


53 


± 





07) 


X 


10^ 


7 


89 


± 





10 





52 


± 





05 


(1 


21 


± 





15) 


X 


109 


325 


± 


19 


109 ± 13 


SDSS 


J230534 


79-010225 


2 


(1 


04 


± 





02) 


X 


10* 


7 


89 


± 





10 





53 


± 





05 


(5 


27 


± 





63) 


X 


10® 


413 


± 


25 


139 ± 13 


SDSS 


J231626 


98+004607 





(5 


13 


± 





02) 


X 


10^ 


7 


90 


± 





10 





52 


± 





05 


(4 


59 


± 





81) 


X 


109 


148 ± 9 


114 ± 8 


SDSS 


J233227 


63-010713 


8 


(1 


20 


± 





02) 


X 


10* 


7 


89 


± 





10 





53 


± 





05 


(3 


58 


± 





45) 


X 


10® 


371 ± 23 


95 ± 9 


SDSS 


J233817 


06-005720 


1 


(9 


74 


± 





12) 


X 


10^ 


7 


89 


± 





10 





53 


± 





05 


(6 


30 


± 





74) 


X 


10® 


424 ± 26 


171 ± 16 


SDSS 


J234110 


13+003259 


5 


(1 


23 


± 





02) 


X 


10* 


7 


89 


± 





10 





53 


± 





05 


(3 


30 


± 





42) 


X 


10® 


334 ± 21 


173 ± 12 


SDSS 


J235138 


85+002716 


9 


(7 


00 


± 





06) 


X 


10^ 


7 


89 


± 





10 





52 


± 





05 


(1 


46 


± 





18) 


X 


109 


319 ± 19 


145 ± 15 



and the resultant distance estimates probably should lie be- 
tween the two extremes presented here. 

The estimated tangential velocity, Vtan, is directly re- 
lated to the distance estimates described above by Kan = 
Dfi, and so it is not surprising that Fig. [TO] exhibits some 
of the same features as Fig. [O] In particular the Bayesian 
values of Kan are correspondingly smaller than those given 
by V07, and only 10 of these 34 sources would satisfy the 
Kan ^ 160 km s~^ criterion as candidate halo WDs if the 
Bayesian estimates were used. Just as before, the Gaussian 
log g prior is probably overly prescriptive, but the essential 
conclusion, that the tangential velocities of these objects 
are lower than estimated by V07, is fairly robust. There are 
also three WDs, discussed further in Section [6. 1.1 1 for which 
the Bayesian estimate of Kan is considerably greater than 
that quoted by V07. They also have the three highest mea- 
sured proper motions of the 34 WDs in the sample, which 
reinforces their identification as halo candidates. Their mea- 
sured proper motions also give some information about their 
direction of motion, which could in principle constrain them 
to have low D, most obviously if this implied an unreason- 
ably high orbital speed along the Solar circle. However this 
is not the case for any of these three sources, and the uncer- 
tainty due to the lack of radial velocity information makes 
it difficult to say more. 

The best way to include this dynamical information 
would be to adopt a velocity prior, in much the same way as 
the distance prior was developed in Section |4]2l As changing 
the distance to a WD with fixed physical properties and a 



fixed spatial velocity, V , would affect both its measured pho- 
tometry and proper motion, it is the only way to correctly 
include the available prior information (i.e., on g, T^s, D and 
V) and thus utili se all the available data simultaneously (cf. 
iHogg et al.ll2005l ). 

All the above results were obtained under the assump- 
tion that the V07 halo candidates have H atmospheres, 
but given that ~ 10 per cent of WDs have He atmo- 
spheres, several of the 34 halo candidates are probably in 
this category. This can be assessed objectively by using the 
model comparison formalism described in Section |3l Com- 
paring En and Ehc reveals that most of the 34 halo can- 
didates are reasonably consistent with either model, with 
I In(iJH) — ln(i?Ho)| ^ 5, although some are decisively bet- 
ter fits to the H models. However, care must be taken in 
interpreting the evidence, as its numerical value itself is of 
little meaning in isolation, especially due to the strong model 
prior that H atmosphere WDs are so more common than He 
WDs. Adopting the model prior tth ~ 0.9 (discussed in Sec- 
tion |4| allows the posterior model probability, Ph, to be cal- 
culated according to Eq. These are shown in Fig. llll and 
compared to the commonly used model comparison statistic 
Ax^ ~ Xmin.Ho " Xmin,H (cf- Sectiou [S]) . As expected, this 
is correlated with Ph, but there is a broad scatter in this 
relation because Ax^ is sensitive only to a single point in 
parameter space (i.e., the best fit) whereas Ph is a properly 
weighted integral over all the possible models. There is also 
a systematic shift due to the strong model prior; if it was 
set to the non-informative value of tth ~ 0.5, all the points 
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Figure 10. The estimated tangential velocity, Vtan, of the 34 
candidate halo WDs as found by V07 and obtained using the 
Bayesian method described here (for H atmosphere models, with 
the Galactic model priors). The horizontal and vertical dotted 
lines show the criterion of Vtan ^ 160 km s~^ used by V07 to 
select their sample of halo candidates. The other two dotted lines 
show the loci Vtan = 0.6 Vtan, grid and Vtan = 0.85 \4an, grid, ex- 
pected for WDs with unconstrained surface gravities for which 
V07 estimate logg = 7.0 and loggr = 7.5, respectively. 



would be higher in Fig. [TT] and the more probable hypothe- 
sis would be simply that with the higher evidence value, as 
seen from Eq. ((Sjl. 

In order to investigate these conclusions further, Table|4] 
gives some of the quantities used in model selection. The 
best fit log g obtained with the uninformative log g prior is 
shown as a proxy for the influence of the informative log g 
prior - the further this is from log g = 7.9, the stronger the 
effect of the g prior will be in the model selection. The best 
fit (for the uninformative prior) is shown as an indicator 
of the infiuence of the data. The evidence values are pre- 
sented to indicate how consistent the data are with each 
model over the entire model space. If neither the g prior 
nor the data definitively picks out an atmosphere model, 
the 10:1 model prior dominates the model selection. Finally, 
the probability of the WD having a H atmosphere is listed, 
along with an indication of the driving force behind this 
conclusion. 

Overall, the implication of these results is that many 
of the candidate halo stars do not actually have the high 
tangential velocities that the data seem to imply. This hy- 
pothesis could be tested by obtaining follow-up observations, 
although higher precision photometry (other than u band) 
would be of little use, as the ambiguities are dominated by 
the fundamental WD model degeneracies, not the quality 
of the existing photometry. Distances based on parallaxes 
would obviously be definitive, but impractical with current 
observational resources. That leaves spectroscopic follow-up, 
although in many cases this is not necessary as 17 of the 34 
halo candidates already have spectroscopic fits presented by 
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''^min.He ''^min.H 

Figure 11. A comparison of the H and He fits to the 34 can- 
didate halo WDs. The difference in Xmin compared to the 
Bayesian probability that the WD has a H atmosphere, Pjj, cal- 
culated assuming the prior probability that any WD has a H 
atmosphere is tth = 0.9. The dashed horizontal line is located at 
Pll = Pho = 0.5; the eight WDs above this line probably have He 
atmospheres, though in no case is this conclusion decisive. Each 
star in the subset of the probable He stars for which spectroscopic 
confirmation exists in the Eisenstein et al. (2006) sample has been 
classified as type DA. 



Table 3. Estimated model parameters of the 17 V07 can- 
didate halo WDs with spectroscopic fits presented by Eisen- 
stein et al. (2006). Three WDs, SDSS J024837.53-003123.9, 
SDSS J024837.53-003123.9 and SDSS J234110.13-I-003259.5, 
were observed twice. 



object type log 9 T„ff 

(K) 



SDSS 


J000557.24-I-001833.2 


DZ 


7.78 ± 0.58 


10110 ± 50 


SDSS 


J001306. 23-1-005506. 4 


DA 


8.13 ± 0.14 


10890 ± 170 


SDSS 


J001838. 54-1-005943. 5 


DA 


7.64 ± 0.13 


21820 ± 760 


SDSS 


J003730. 58-001657.8 


DA 


7.93 ± 0.38 


9130 ± 160 


SDSS 


J003813. 53-000128. 8 


DA 


8.03 ± 0.11 


12140 ± 400 


SDSS 


J005906. 78-1-001725. 2 


DZ 


8.53 ± 0.38 


11180 ± 270 


SDSS 


J010207.24-003259.7 


DA 


8.24 ± 0.08 


11050 ± 100 


SDSS 


J010225. 13-005458.4 


DA 


8.95 ± 0.06 


18790 ± 830 


SDSS 


J014247.10-I- 005228. 4 


DA 


8.82 ± 0.11 


14040 ± 450 


SDSS 


J024837.53-003123.9 


DA 


7.81 ± 0.15 


9080 ± SO 






DA 


8.11 ± 0.17 


9500 ± 110 


SDSS 


J025325. 83-002751. 5 


DA 


7.70 ± 0.03 


19200 ± 200 


SDSS 


J025531. 00-005552. 8 


DA 


7.72 ± 0.20 


12720 ± 640 






DA 


7.74 ± 0.15 


12830 ± 440 


SDSS 


.7211928.44-002632. 9 


DA 


7.87 ± 0.10 


15660 ± 410 


SDSS 


J215138. 09+003222. 3 


DA 


8.23 ± 0.45 


7820 ±210 


SDSS 


J233227.63-010713.8 


DA 


7.60 ± 0.20 


14420 ± 1320 


SDSS 


J233817.06-005720.1 


DA 


8.07 ± 0.33 


10380 ± 280 


SDSS 


J234110.13-)-003259.5 


DA 


7.90 ± 0.08 


13380 ± 330 






DA 


7.80 ± 0.09 


12830 ± 280 



lEisenstein et al.l ()2006h . These spectroscopic fits (only log 
and Teff, but not D) are given in Table |31 and compared to 
the Bayesian estimates of Teff in Fig. 1121 

The surface gravity estimates obtained with the non- 
informative prior reveal that the photometric data for the 
halo candidates are most consistent with low g, as expected 
given the combination of solid angle degeneracy and high- 
Kan selection. There is, however, a suggestion that the true 
surface gravities are systematically higher than the fiducial 
value of log g ~ 7.9, a notion that is supported by the fact 
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Table 4. Model selection quantities for the 34 V07 candidate halo WDs, where columns 2 and 3 are the best-fit logg values 
estimated using the uninformative g prior. UKIDSS data have been used where available. The final column states whether the 
model selection result is driven by the data, the informative g prior, or the 10:1 model prior. 



object logSH logSHc Xmin.H Xmin.Ho -^H notcs 



SDSS 


J000244. 03+010945. 8 


7 


57 


7.83 


16.5 


14.1 


1 


2 


X 


j(,-02 


1 


2 


X lO-"! 





47 


data, and g prior 


SDSS 


J000557. 24+001833. 2 


7 


86 


7.89 


2.8 


4.1 


5 


6 


X 




2 


3 


X 10-02 





96 


data and model prior 


SDSS 


J001306. 23+005506. 4 


7 


55 


7.84 


14.8 


12.6 


3 


1 


X 




4 


7 


X 10-02 





37 


data and g prior 


SDSS 


J001518. 33+010549. 2 


7 


38 


7.94 


13.3 


8.9 


1 


8 


X 


10-02 






1.3 





11 


data and g prior 


SDSS 


J001S3S. 54+005943. 5 


7 


57 


7.97 


11.1 


10.3 


7 


1 


X 


io-°i 






4.5 





58 


model prior 


SDSS 


J002951. 64+005623. 6 


7 


74 


7.87 


10.5 


8.6 


7 


1 


X 


10-01 






4.4 





59 


model prior 


SDSS 


J003054. 06+001115. 6 


7 


90 


7.95 


23.3 


30.3 


6 





X 


10-04 


1 


3 


X 10~05 


1 


00 


data, neither model good fit 


SDSS 


J003730. 58— 001657.8 


7 


81 


7.88 


0.9 


3.4 


1 


9 


X 


10+01 






4.5 





97 


data and model prior 


SDSS 


J003813. 53— 000128.8 


7 


61 


7.77 


3.4 


3.9 


3 


3 


X 


io-°i 




5 


X 10-"! 





84 


m.odel prior 


SDSS 


J004214. 88+001 135. 7 


7 


68 


7.79 


12.4 


11.4 


2 


1 


X 


io-«i 


9 


7 


X 10-01 





66 


model prior 


SDSS 


J005906. 78+001725.2 


7 


73 


7.84 


6.0 


4.8 


3 


9 


X 


io-«i 






2.1 





62 


m.odel prior 


SDSS 


J010129. 81— 003041.7 


7 


54 


7.86 


16.0 


15.2 


3 


2 


X 


10-02 


1 


8 


X 10-"! 





62 


m.odel prior 


SDSS 


J010207. 24— 003259.7 


7 


93 


7.79 


0.5 


1.1 






1 


.1 


9 


6 


X W-°^ 





91 


data, g and model prior 


SDSS 


J010225. 13—005458.4 


7 


65 


7.96 


0.9 


0.2 


3 


9 


X 


10+01 


1 


2 


X 10 + 02 





75 


model prior 


SDSS 


J014247.10+005228.4 


7 


46 


7.87 


2.1 


0.4 






6.6 


5 


7 


X 10 + 01 





51 


model prior 


SDSS 


J015227. 57-002421.1 


7 


81 


7.79 


1.5 


1.7 


1 


8 


X 


10+01 


1 


6 


X lO+t'l 





91 


model prior 


SDSS 


J020241. 81-005743.0 


7 


82 


7.89 


9.3 


11.2 


4 


6 


X 


10-02 


1 


2 


X 10-02 





97 


data and model prior 


SDSS 


J020729. 85+000637. 6 


7 


78 


8.61 


8.4 


23.4 


6 


6 


X 


10-03 


8 


6 


X 10-"** 


1 


00 


data 


SDSS 


J024529. 69-004229. 8 


7 


43 


7.86 


8.7 


6.0 


7 


2 


X 


10-01 


1 


2 


X 10 + 01 





35 


data and g prior 


SDSS 


J024837. 53-003123. 9 


8 


08 


8.10 


15.3 


31.6 


4 


2 


X 


10-05 


2 


6 


X 10-"** 


1 


00 


data 


SDSS 


J025325. 83-002751. 5 


7 


13 


7.93 


97.3 


118.3 


4 


1 


X 


10-2" 


6 


2 


X 10-26 


1 


00 


data, neither model is a good fit 


SDSS 


J025531. 00-005552. 8 


7 


64 


7.91 


5.7 


3.7 






2 


.2 


1 


6 


X 10 + 01 





55 


model prior 


SDSS 


J030433. 61-002733. 2 


7 


74 


8.22 


15.6 


31.0 


7 


1 


X 


10-03 


1 


3 


X 10-06 


1 


00 


data and g prior 


SDSS 


J211928. 44-002632. 9 


7 


41 


8.06 


7.7 


4.7 


7 


8 


X 


10-01 


1 


3 


X 10 + 01 





36 


data and g prior 


SDSS 


J213641. 39+010504. 9 


7 


54 


7.92 


13.9 


11.0 


4 





X 


10-02 


5 


8 


X 10-01 





38 


data and g prior 


SDSS 


J215138. 09+003222. 3 


7 


99 


8.09 


7.1 


22.6 


2 


3 


X 


10-01 


8 





X 10 -OS 


1 


00 


data, g and model prior 


SDSS 


J223808. 18+003247. 9 


7 


64 


8.43 


40.2 


75.4 


1 


7 


X 


10-09 


9 


2 


X 10-1** 


1 


00 


data, neither model is a good fit 


SDSS 


J223815. 97-011336. 9 


7 


55 


7.77 


17.4 


10.9 


4 


2 


X 


10-03 


1 


8 


X 10-"1 





17 


data and g prior 


SDSS 


J230534. 79-010225. 2 


7 


71 


7.81 


9.0 


7.9 


7 


6 


X 


io-"i 






2.4 





74 


model prior 


SDSS 


J231626. 98+004607.0 


7 


88 


8.34 


7.4 


0.6 


1 


4 


X 


10-02 


1 


1 


X 10-01 





55 


g and model prior 


SDSS 


J233227. 63-010713. 8 


7 


64 


7.84 


4.4 


2.1 






1 


.7 


1 


1 


X 10 + 01 





57 


model prior 


SDSS 


J233817. 06-005720.1 


7 


94 


7.85 


17.5 


21.4 


7 


5 


X 


10-03 


1 


9 


X 10-"^ 





97 


data 


SDSS 


J234110. 13+003259. 5 


7 


54 


7.86 


3.5 


1.3 






1 


.8 


2 


4 


X 10 + 01 





40 


data and g prior 


SDSS 


J235138. 85+002716. 9 


7 


67 


7.94 


10.5 


16.2 


1 


7 


X 


io-"i 


8 


4 


X 10-"3 





99 


data and model prior 



that the available spectroscopic log g fits are similarly dis- 
crepant. The simplest explanation for this is that the Gaus- 
sian log g prior given in Eq. (|10p is overly prescriptive. Cor- 
recting the tails of the prior would not alter the best fit 
values, but would proportionally increase the tails of log g 
posteriors. Given that the solid angle degeneracy ensures 
that a broad range of log g values are consistent with the 
data, the photometric measurements add minimal new in- 
formation, and so the posterior peaks at the prior value for 
each WD in turn. Some particular cases of this phenomenon 
are discussed in Section [6. 1.11 



The spectroscopic and photometric temperature esti- 
mates shown in Fig. [12] are reasonably consistent with each 
other, but the Bayesian values are systematically lower than 
those from the spectroscopic fits. This is not the result of 
the steep Tcs prior given in Eq. (|lip . but is another, as 
yet unidentified, systematic difference between the two ap- 
proaches. The natural assumption would be that the spec- 
troscopic fits are correct given the far greater available in- 
formation; however, as can be seen in Table [S] many of 
these WD have temperatures of T^s < 11, OOP K for whic h 
the spectroscopic models are biased ([Eisenstein et al.]|2006t ). 
A final possibility is, once again, a statistical explanation: 
whereas there is no reason that the spectroscopic data would 
be at the margins of the noise distribution, it is plausible 
that the photometric data on the same objects are biased, 
as it was these data that were used to select the sample. 



(b) 



5000 



10* 



2x10* 
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eff ,spectroscop: 



.e (K) 



Figure 12. The estimated effective temperature, Tj,g, of the 17 
of the 34 candidate halo WDs identified by V07 for which there 
are spectroscopic fits in Eisenstein et al. (2006). Here the values 
estimated from spectra are compared to the Bayesian estimates 
obtained from the H atmosphere models using the informative 
Galactic model prior. 
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Figure 13. Marginalised posterior distributions of the model parameters for the candidate halo WD SDSS J010207. 24— 003259.7. The 
point with the error bar in the left panel is the Eisenstein et al. (2006) fit to the SDSS spectrum of this WD. For both the non-informative 
prior (grey) and the informative prior (orange) the coloured regions enclose 68 per cent and 95 per cent of the posterior probability in 
these two parameters. The open triangle in the right panel is the Harris et al. (2006) fit assuming logg = 8.0. 



6.1.1 Notes on individual objects 

The results for most of the 34 candidate halo WDs are ex- 
plained by the general arguments given above, but it is useful 
to investigate a few of sources in greater detail. 

Of particular interest are the three objects for 
which spectrosc o pic s urface gravity estimates given by 
lEisenstein et al.l l|2006l ) are significantly higher than the 
fiducial value of logs - 7.9: SDSS J005906.78-h001725.2, 
with logp = 8.53 ±0.38; SDSS J010225.13-005458.4, with 
logg = 8.95 ±0.06; and SDSS J014247.10±005228.4, with 
logg — 8.82 ± 0. 11. The first of these was identified as a 
DZ type WD by iKleinman et al] (|2004l ) . and it is perhaps 
not surprising that the DA and DB atmosphere fits pre- 
sented here and in V07 do not match the spectroscopic es- 
timate. The latter two objects were identified as DA type 
WDs bv lEisenstein et al.1 (|2006l ). There are no UKIDSS mea- 
surements for these two stars, and the goodness-of-fit to the 
SDSS photometry shows nothing unusual. The implication is 
that these objects have intrinsically unusual surface gravity, 
and that the photometry is insufficient to add any informa- 
tion to the expectations of the log g prior. It further follows 
that, because of the direction of the degeneracy between g 
and D, a higher value of log g would imply a significantly 
lower distance, and hence, lower Kan- 

There are also three WDs seen in Fig. [10] (b) with strik- 
ingly high Bayesian estimates for Kan compared with the 
V07 estimates: SDSS J010207.24-003259.7, with Kan ^ 
334 kms-\ SDSS J024837.53-003123.9, with Kan ^ 
328 kms-\ and SDSS J223808.18±003247.9, with Kan ^ 
289 km s"^ There are UKIDSS observations for all three 
st ars and, fortuitously, sp ectroscopic log g estimates given 
bv lEisenstein et al] (|2006l ) for the first two, the marginalised 
posterior distributions for which are shown in Figs. [131 and 
[M respectively For SDSS J010207.24-003259.7 the spec- 
troscopic estimate of logg = 8.24 ± 0.08 implies that it 
is probably not a viable halo candidate. The spectroscopic 
log estimates for SDSS J024837.53-003123.9 are consistent 
with the Bayesian estimates, and hence imply this WD may 
have genuinely unusual dynamics. The fit for the third ob- 
ject, SDSS J223808.18±003247.9, which does not have spec- 



J024837.53-003123.9 




0.85 0.90 0.95 



1.00 



Figure 14. Same as Fig. I13f a^. but for the candidate halo WD 
SDSS J024837.53-003123.9. The two points with errors are the 
Eisenstein et al. (2006) fits to the two SDSS spectra of this WD. 



troscopic data, is disastrously bad for both H and He models, 
with or without UKIDSS data. 

S PSS J024837.53- ( )0312 3.9 has duplicate observations 
in the [Eisenstein et al.l ([2006) sample (see the left panel of 
Fig. I14p . and it is interesting to note that the error bars 
given for Toff for these repeated observations are marginally 
inconsistent with each other at face value (as well as being 
systematically higher than the Bayesian estimate, as men- 
tioned above). The temperature uncertainties from the pho- 
tometry are comparably precise to those from spectroscopy, 
although the surface gravity is only constrained by the spec- 
troscopic data. 

Although the model constraints on 
SDSS J030433.61-002733.2 are not interesting in them- 
selves it is important to note that the UKIDSS data 
associated with this object in V07 is results from an 
incorrect cross- match to a much cooler star 3'.' 5 away. 
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Table 5. Estimated H atmosphere model parameters of the 24 V07 candidate ultra-cool WDs, obtained using the 
informative prior defined in Eq. l[9]l. 



object T^ff log g M T D Vta„ 

(K) (Mq) (yr) (pc) (kms-l) 



SDSS 


J001107 


57—003102 


8 


(3 


64 


± 





04) 


X 


10^ 


7 


36 


± 





06 





27 


± 





02 


(4 


28 ± 0.37) 


X 


10^ 


105 ± 3 


28 ± 3 


SDSS 


J004S43 


28-003820 





(3 


72 


± 





05) 


X 


10^ 


7 


53 


± 





09 





34 


± 





04 


(5 


07 ± 0.57) 


X 


10^ 


146 ± 6 


37 ± 6 


SDSS 


J012102 


99 — 003833 


6 


{3 


93 







04) 


X 


10^ 


7 


75 







07 





45 







03 


(6 


24 ± 0.66) 


X 


10^ 


67 ± 3 


39 ± 2 


SDSS 


J013302 


17+010201 


3 


(3 


62 


+ 





10) 




10^ 


7 


62 


+ 





10 





38 


+ 





05 


(6 


13 ± 0.75) 




10^ 


145 ± 7 


34 ± 6 


SDSS 


J030144 


09-004439 


5 


(3 


62 


± 





05) 


X 


10^ 


7 


77 


± 





05 





46 


± 





03 


(7 


28 ± 0.52) 


X 


109 


50 ± 2 


132 ± 4 


SDSS 


J204332 


97+011436 


2 


(3 


56 


± 





07) 


X 


10^ 


7 


62 


± 





09 





38 


± 





05 


(6 


35 ± 0.73) 


X 


109 


139 ± 7 


38 ± 9 


SDSS 


J205010 


17+003233 


7 


(3 


85 


± 





14) 


X 


10^ 


7 


85 


± 





10 





50 


± 





05 


(7 


55 ± 1.14) 


X 


109 


92 ± 7 


24 ± 6 


SDSS 


J205132 


05+000353 


6 


(3 


46 


± 





04) 


X 


10^ 


7 


05 


± 





04 





17 


± 





01 


(3 


18 ± 0.17) 


X 


10^ 


69 ± 2 


331 ± 9 


SDSS 


J210742 


26-002354 


1 


(3 


49 


± 





08) 


X 


10^ 


7 


77 


± 





09 





46 


± 





05 


(7 


86 ± 0.81) 


X 


10^ 


109 ± 5 


42 ± 6 


SDSS 


J212216 


01-010715 


2 


(4 


09 


± 





06) 


X 


10^ 


7 


91 


± 





10 





53 


± 





05 


(7 


53 ± 1.23) 


X 


10^ 


105 ± 7 


84 ± 7 


SDSS 


J212930 


25-003411 


5 


(3 


93 


± 





06) 


X 


10^ 


7 


98 


± 





10 





57 


± 





06 


(8 


56 ± 0.98) 


X 


109 


76 ± 6 


154 ± 11 


SDSS 


J214108 


42+002629 


5 


{3 


60 


± 





10) 


X 


10^ 


7 


60 


± 





10 





37 


± 





05 


(6 


03 ± 0.73) 


X 


10^ 


150 ± 8 


38 ± 8 


SDSS 


J220455 


03-001750 


6 


(3 


76 


± 





14) 


X 


10^ 


7 


97 


± 





10 





57 


± 





06 


(8 


97 ± 1.06) 


X 


10^' 


81 ± 7 


41 ± 4 


SDSS 


J223105 


29+004941 


9 


(3 


72 


± 





04) 


X 


10^ 


7 


49 


± 





07 





32 


± 





03 


(4 


78 ± 0.41) 


X 




113 ± 3 


61 ± 3 


SDSS 


J223520 


19-003623 


6 


(3 


74 


± 





10) 


X 


10^ 


7 


91 


± 





09 





53 


± 





05 


(8 


41 ± 1.17) 


X 


10^ 


110 ± 7 


73 ± 6 


SDSS 


J223715 


32-002939 


2 


(3 


30 


± 





06) 


X 


10^ 


7 


72 


± 





08 





43 


± 





04 


(8 


05 ± 0.60) 


X 


10^ 


99 ± 4 


30 ± 3 


SDSS 


J223954 


07+001849 


2 


(4 


22 


± 





06) 


X 


10^ 


7 


79 


± 





09 





47 


± 





05 


(6 


00 ± 0.91) 


X 


10^ 


79 ± 4 


41 ± 3 


SDSS 


J224206 


19+004822 


7 


(3 


40 


± 





03) 


X 


10^ 


8 


29 


± 





06 





77 


± 





04 


(1.07 ± 0.01) 


X 


IQlO 


24 ± 1 


18 ± 1 


SDSS 


J224845 


93-005407 





(4 


02 


± 





05) 


X 


10^ 


7 


93 


± 





10 





54 


± 





06 


(7.90 ± 1.22) 


X 


10^ 


91 ± 6 


88 ± 7 


SDSS 


J225244 


51+000918 


6 


(3 


73 


± 





14) 


X 


10^ 


7 


94 


± 





10 





55 


± 





06 


(8 


70 ± 1.25) 


X 


109 


96 ± 9 


39 ± 5 


SDSS 


J232115 


67+010223 


8 


(4 


46 


± 





03) 


X 


10^ 


7 


94 


± 





10 





55 


± 





06 


(6 


79 ± 1.07) 


X 


10^ 


53 ± 3 


68 ± 5 


SDSS 


J233055 


19+002852 


2 


(4 


23 


± 





03) 


X 


10^ 


7 


89 


± 





08 





52 


± 





04 


(7.08 ± 0.96) 


X 


109 


49 ± 2 


39 ± 2 


SDSS 


J233818 


56-004146 


2 


(3 


70 


± 





06) 


X 


10^ 


7 


72 


± 





09 





43 


± 





04 


(6 


66 ± 0.82) 


X 


109 


122 ± 6 


82 ± 6 


SDSS 


J234646 


06-003527 


6 


(3 


70 


± 





06) 


X 


10^ 


7 


55 


± 





09 





35 


± 





04 


(5.24 ± 0.60) 


X 


109 


143 ± 6 


30 ± 5 



Table 6. Estimated He atmosphere model parameters of the 24 V07 candidate ultra-cool WDs, obtained using the 
informative prior defined in Eq. JgJ. 



object T„ff log g M T D Vta„ 

(K) (Mq) (yr) (pc) (kms-1) 



SDSS 


J001107 


57-003102 


8 


(4 


49 


± 





04) 


X 


10^ 


7 


99 


± 





11 





57 


± 





06 


(6 


68 


± 





78) 


X 


109 


106 ± 9 


29 ± 4 


SDSS 


J004843 


28-003820 





(4 


42 


± 





04) 


X 


10^ 


7 


93 


± 





10 





53 


± 





05 


(6 


41 


± 


1 


00) 


X 


109 


154 ± 11 


39 ± 7 


SDSS 


J012102 


99-003833 


6 


(4 


65 


± 





02) 


X 


10^ 


7 


92 


± 





10 





53 


± 





05 


(5 


83 


± 


1 


00) 


X 


109 


83 ± 5 


48 ± 4 


SDSS 


J013302 


17+010201 


3 


(4 


35 


± 





05) 


X 


10^ 


7 


91 


± 





10 





53 


± 





05 


(6 


42 


± 


1 


07) 


X 


109 


159 ± 12 


38 ± 7 


SDSS 


J030144 


09-004439 


5 


(4 


60 


± 





02) 


X 


10^ 


7 


90 


± 





10 





52 


± 





05 


(5 


74 


± 


1 


05) 


X 


109 


70 ± 4 


185 ± 12 


SDSS 


J204332 


97+011436 


2 


(4 


33 


± 





05) 


X 


10^ 


7 


90 


± 





10 





52 


± 





05 


(6 


39 


± 


1 


09) 


X 


109 


155 ± 12 


43 ± 10 


SDSS 


J205010 


17+003233 


7 


(4 


60 


± 





06) 


X 
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6.2 Candidate ultra-cool white dwarfs 

The 24 candidate ultra-cool WDs identified by V07 were 
analysed using tiie teciiniques described in Section [2] and 
Section [31 giving parameter estimates and evidence values 
for both H and He atmosphere models. The results obtained 
using the informative prior defined in Section|4]are presented 
in Tables[S]and[51 The Bayesian parameter fits are compared 
to those obtained by V07 in Fig.[T5]for the informative prior. 

For the ultra-cool WDs, the temperature comparison 
shown in Fig. [15] is the critical test as it was on the basis of 
Tcff.grid that these 24 objects were selected. The discretisa- 



tion on the horizontal axis is again due to the grid-based fit- 
ting of V07. The slight tendency for the Bayesian fits to have 
a higher Tcff is due to the selection effect that resulted in 
the more pronounced distance bias discussed in Section [6. II 
The various distance estimates for the ultra-cool WDs 
are also in broad agreement. Comparing the results obtained 
using the non-informative prior with those provided by V07 
reveals only a small scatter due to the best-fit models lying 
between the model grid points. The Bayesian distance esti- 
mates obtained from the informative prior and those given 
by V07 also exhibit systematic differences, with the former 
tending towards the average value of D ~ 100 pc. 
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Figure 15. The estimated effective temperature, Teffi of the 24 
candidate ultra-cool WDs as found by V07 and obtained using the 
Bayesian method described here (for H atmosphere models, with 
the Galactic model priors). The selection criterion as an ultra-cool 
WD, Teff ^ 4000 K, adopted by V07 is shown as the vertical and 
horizontal dotted lines. One source, SDSS J232115.67-f 010223.8, 
is listed in this sample despite an estimated effective temperature 
above this nominal limit. 

Given that the Bayesian distance estimates are in agree- 
ment with those of V07, it is not surprising that the same 
holds for the tangential velocity estimates. The only signif- 
icant exception to this is SDSS J205132.05-f 000353.6, for 
which the H fit shown is probably inappropriate, confirming 
the findings of V07 that this WD has a He atmosphere. This 
object is discussed further in Section [6.2.11 

The model comparison results for the ultra-cool sample 
are summarised in Table [T] The most striking result is that 
18 out of the 24 candidate ultra-cool WDs have a high prob- 
ability of being He stars, despite the strong model prior that 
H stars should be ~ 10 times more common. The same con- 
clusion is also implied by the optical data alone, although 
it is perhaps more revealing that several stars are poor fits 
to both the H and He models. It is also interesting to note 
that the best-fit temperatures of the He stars are typically 
higher than for H in this region of parameter space. Thus it 
is possible that selecting WDs on the basis of having a low 
Tcff from H fits results in also selecting He stars with less 
extreme temperatures. 



6.2.1 Notes on individual objects 

Three stars in the sample of candidate ultra-cool WDs 
have been spectroscopically confirmed as being of type DC. 
SDSS J224206+004822 and SDSS J233055-H002852 were 
initially discovered by iKilic et"aL I (12006 ) who, from pho- 
tometric fits with log g = 8.0, found TaS — 3400 K and 
Tefs = 4100 K, respectively. Applying the Bayesian ap- 
proach described in Section 12.21 and using the informative 
log g prior given in Section |4] gives Tag = 3400 ± 30 K for 
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Figure 16. Photometry in the g, r, i and 2bands (solid dots) 
of the probable He WD SDSS J205132-I-000353, compared to the 
best-fit H models (open squares) and He models (open triangles) 
in the n, g, r, i, z, J, H and K bands. The photometric errors 
are ~ 0.02, and smaller than the symbol size. 



SDSS J224206-h004822 (where H is a better fit but nei- 
ther model fits the data well), and f^g = 4760 ± 20 K for 
SDSS J233055-f 002852 (where He is clearly a better fit) 
SDSS J232115.67-H010223.8 was reported by ICaroUo et all 
(|2006[ ) as being of type DC; the Bayesian results imply 
that this star most probably has a He atmosphere, with 
Teff = 4870 ± 20 K. 

Only one of the 24 ultra-cool WD candidates identified 
by V07 was found to definitely have a He-dominated at- 
mosphere. Given the strong degeneracies in the photomet- 
ric WD fitting, it was possible that this conclusion came 
about only due to the sparse grid of H models 'strad- 
dling' the good fit locus while one of the He models fortu- 
itously landed on it. Nonetheless, these inferences are sup- 
ported by the more thorough MCMC exploration of the 
model space, although the fact that neither atmospheric 
species results in a good fit (as shown in Fig. I16p leaves 
the possibility that SDSS J205132-I-000353 is not a WD at 
all. This possibility has been corroborated by Kilic (pri- 
vate communication), who did not detect an appreciable 
proper motion from multi-epoch astrometry of this source, 
implying it is more likely to be a quasar, or other extra- 
Galactic object. These observations are at odds with proper 
motion of /i = (1.02 ± 0.01) arcsec yr~^ quoted by V07; 
however there are two other stars within a few arcsec of 
SDSS J205132+000353, and so it is possible t h at th e au- 
tomated matching algorithm of iBramich et af] (|2008| ) was 
unable to process this complicated case. 

SDSS J205132+000353 and SDSS J212930.25-003411.5 
have very high Kan values of 262 ± 21 km s'^ and 224 ± 
15 km s~^, respectively. The fact that the former is not a 
good fit to either model atmosphere probably invalidates 
the estimate of Kan, but the latter remains a good halo 
candidate. 
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Table 7. Model selection quantities for the 24 V07 candidate ultra-cool WDs, where columns 2 and 3 are the best fit logg values 
estimated using the uninformative g prior. UKIDSS data have been used where available. The final column states whether the 
model selection probability of H versus He is driven by the data, the informative g prior, or the 10:1 model prior. 
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7 CONCLUSIONS 

WDs can be identified reliably using just photometric and 
astrometric information, but the utility of large samples 
generated in this way is limited if WDs cannot also be 
characterised without recourse to spectroscopy. The degree 
to which this is possible has been assessed here by using 
Bayesian methods to include vital prior information in the 
parameter estimation process. For most WDs observed with 
SDSS precision photometry (or per cent-level photometry, 
but only in bands redder than u) only their effective tem- 
perature and solid angle can be constrained unambigously. 
However this very generic degeneracy can be broken by util- 
ising the prior information (from spectroscopic modelling) 
that most WDs have surface gravities of 7.7 < log g < 8.1. 

The power of this approach was demonstrated by ap- 
plying it to the samples of candidate halo and ultra-cool 
WDs generated by V07 from muhi-epoch SDSS Stripe 82 
data. This confirms their identification of ultra-cool stars 
(as the photometric data constrain the temperature well) 
but implies that most of the candidate halo stars are mem- 
bers of the more numerous thick disk population with un- 
usual observed colours caused by photometric noise. This is 
a particular case of a very general result, that it is almost 
impossible to reliably identify a small number of outliers if 
there are sufficiently many members of the more dominant 
population to scatter to the parameter space of interest. 

The one situation in which photometric data are likely 
to give definitive WD parameter estimates without prior in- 
formation is if per cent-level photometry is available in the 
u, g and r (or similar) bands for H at mosphere WDs wit h 
10000 K < Tcff < 18000 K. As shown bv llvezic et aD l|2007h . 
the loci of observed colours of WDs in SDSS Stripe 82 ex- 
hibit an obvious separation between H and He models (see 
Fig. [T]), and the next stage of this work will be to analyse a 
sample of such sources. This should also have the additional 
benefit of giving a more accurate surface gravity prior than 



used here (in which the numbers of WDs with non-fiducial 
surface gravities are almost certainly under-estimated). An- 
other important extension to this work will be to add a 
velocity prior from Galactic dynamics; for most WDs this 
will not be important, but it will be critical to obtaining 
the most stringent constraints on potential halo objects. An 
added side-benefit of being able to identify halo objects dy- 
namically will be to place empirical constrain ts on the tem- 
perature and age distribution of halo WDs (cf. lHanse nl l200ll : 
iPauh et al.ll2006l '). Further, while the effect of the ISM has 
not been included in this work (in part to follow the as- 
sumptions of the V07 analysis) , it would be straightforward 
to include a full ISM model or Galactic extinction template 
as described in Section (2] 

These results have particularly relevance to up-coming 
large -scale deep optical surveys such as the SkyMapper sur- 
vev llKeller et al.l I2OO7I '). the Dark Energy Survey (DBS; 



iFlaugher et al. 20061) . the Panoramic Surv ey Telescope 
And Rapid Response System (Pan-STARRS; iKaiser et al.l 
I200I) and the Large Synoptic Survey Telescope (LSST; 
Ivezic et al. ,200^ ). In the context of WD studies, a critical 
distinction between these projects is that only SkyMapper 
and LSST plan to observe in the u band. With LSST ex- 
pected to reac iS u ~ 26 over half the sky, it should yield 
a sample of ~ 10^ WDs with sufficiently precise photomet- 
ric and astrometric measurements to give accurate distances 
and tangential velocities. Combined with the methods pre- 
sented here, these data will provide an unprecedented un- 
derstanding of the Galaxy's WD population. 



Typical LSST exposures will detect point-sources with 
a signal— to-noise ratio of 5 to u ~ 24; the depth 
quoted in the main text comes from co-adding the 
^ 60 scans LSST will make during its lifetime. See 
http://www.lsst.org/Science/lsst_baseline.shtml for more 
details. 
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